Tight-binding model for iron pnictides 
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We propose a five-band tight-binding model for the Fe-As layers of iron pnictides with the hopping 
amplitudes calculated within the Slater-Koster framework. The band structure found in DFT, 
including the orbital content of the bands, is well reproduced using only four fitting parameters 
to determine all the hopping amplitudes. The model allows to study the changes in the electronic 
structure caused by a modification of the angle a formed by the Fe-As bonds and the Fe-plane and 
recovers the phenomenology previously discussed in the literature. We also find that changes in a 
modify the shape and orbital content of the Fermi surface sheets. 



I. INTRODUCTION 

Since the discovery of high temperature superconduc- 
tivity in iron pnictides 1 ^ a lot of attention has been de- 
voted to their understanding. Iron pnictides are layered 
materials with arsenic (or another pnictogen) atoms at 
the center of the Fe plaquettes, out of plane and arranged 
upwards and downwards in a checkerboard form (see 
Fig. [J) in tetrahedral configuration. Fe-As bonds form an 
angle a with the Fe-plane, called in the following Fe-As or 
iron-pnictogen angle, which differs among compounds 3,4 
and depends on doping 3 -!^ or applied pressure. 7 A pos- 
sible connection between the value of a, the critical tem- 
perature and electronic properties has been discussed by 
several authors^^ Recently Kuroki et aii°- have pro- 
posed that the pnictogen height above the Fe plane is 
the key factor that determines both T c and the form of 
the superconducting gap. 

From a DFT point of view iron superconductors have 
multiband character mostly due to Fe d-orbitalsj 11 ' 12 ' 13 
In the iron (unfolded) Brillouin zone 1 ^ the Fermi surface 
consists of electron pockets at the X and Y points, two 
hole pockets in T and a hole pocket at M, in reason- 
able agreement with de Haas van Alphen 1 ^. experiments 
in the non-magnetic state. Angle-resolved photoemis- 
sion (ARPES) measurements give also evidence of Fermi 
pockets at these symmetry pointsi 16 ' 17 ' 18 ' 19 ' 20 Interband 
scattering between electron and hole pockets has been 
proposed as a mechanism for superconductivity 1 ' 21 In 
this context the importance of nesting for superconduc- 
tivity and magnetism is discusse d 22 ' 23 . More recently, 
the relevance of the anisotropic orbital weight of each 
Fermi pocket in determining the symmetry of the super- 
conducting order parameter has been emphasize d 10 ' 24 ' 25 . 
Electron pockets at X and Y have respectively yz/xy 
and zx/xy origin while the hole pockets in T arise from 
zx and yz orbitals 2 ^. Due to closeness of two hole bands 
and different dependence of their energy on a the or- 
bital character of the pocket in M switches between xy 
or 3z 2 



r ' r '> \ X (M') 
kx / 



FIG. 1: Top figure: Sketch of the lattice structure. Fe- 
As bonds form an angle a with the Fe-plane which changes 
among compounds, with doping and with pressure. Bottom 
figure: On the left, in a top view of the Fe-As layer, the real 
(extended) and the Fe unit cells are shown in dashed and solid 
lines respectively. The X and Y axis of the Fe unit cell, used 
in the paper and, shown with arrows, are directed along the 
Fe bonds. On the right the experimental (folded) Brillouin 
zone is shown with dashed lines. Its symmetry points are 
denoted with primed letters. The extended Brillouin zone, 
used in the paper, is delimited by solid lines. It is double- 
sized and rotated 45 degrees with respect to the experimental 
Brillouine zone. Bands and Fermi pockets at F and M in the 
extended Brillouin zone and discussed through the text will 
appear experimentally at V . 



of any theoretical treatment in a lattice. Initial attemps 
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or 



" depending on the value of a . ' Experimen- 
tally, the orbital content can be studied by changing the 
polarization of the light used in ARPES.22 1 22i2L22 

A good tight-binding model is the basic building-block 



tried to describe the iron pnictides using two 
three-orbital^ 5 , models. Several proposals based on the 
symmetry properties of zx and yz (and xy as third or- 
bital) were put forward to describe the bands close to the 
Fermi level. However crystal field splittings among the Fe 
d-orbitals are small compared with the bandwidth result- 
ing in strong hybridization of all d-orbitals. At present it 
is believed that inclusion of all five Fe d-orbitals is nec- 
essary to obtain a good description of the properties of 
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iron pnictide o 36 i 37 ' 38 . The placement of As at the center 
of the plaquettes suggests that hopping between Fe atoms 
to second nearest neighbors cannot be disregarded. 

In this paper we propose a five-orbital tight bind- 
ing model to describe the Fe-As layers with the hop- 
ping amplitudes calculated within the Slater-Koster 
framework. 39 Compared to DFT tight-binding fits, the 
procedure presented here greatly reduces the number of 
fitting parameters necessary to calculate the bands and 
allows to study changes in the iron-pnictogen angle a. 
We show that the bands close to the Fermi level can be 
described giving all the hopping amplitudes in terms of 
just four parameters. The agreement between our results 
and DFT predictions extends to the orbital weight of each 
band. We also reproduce the switch in orbital character 
of the hole pocket in M, when a varies. Furthermore, we 
predict that changes in a can induce modifications in the 
shape and orbital content of the Fermi pockets, includ- 
ing the disappearance of the hole pockets in T when the 
tetrahedron is elongated. Within the present theoreti- 
cal understanding these results have strong implications 



in the superconducting and magnetic properties of these 
compounds i 10 ' 24 i 25 



II. TIGHT-BINDING MODEL 

We construct a tight-binding model to describe the 
band structure of the FeAs layers including the five Fe 
d-orbitals. Arsenic atoms only enter in the model indi- 
rectly via the Fe-Fe hopping amplitudes. Indirect hop- 
ping via arsenic is treated to second order in perturba- 
tion theory Direct hopping between Fe atoms is also 
included. Hopping is restricted to first and second near- 
est Fe neighbors. Both the mathematical form of the 
Hamiltonian and the hopping amplitudes are computed 
within the Slater-Koster formalism^ 9 -. We take X and Y 
directed along the Fe-bonds (see Fig[T]). The same axis 
convention applies for the orbitals, i.e. x 2 — y 2 orbital 
lobes are directed along the Fe-Fe bonds. Under these 
assumptions the Hamiltonian is given by: 
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Here m,n refer to lattice sites, 7,/? are the orbital in- 
dices, a the spin and \i the chemical potential. Only 
the first sum runs through all the orbitals. Brackets and 
parentheses restrict the orbitals to which the other sums 
apply. In particular, (7 7^ f3) is restricted to the pairs 
{li P} = {yz,zx}, {xy,3z 2 - r 2 }, [7 ^ (3} to the pair 
{7,/?} = {3z 2 -r 2 ,x 2 -y 2 }, (7^/3) to the pairs {7,/?} = 
{yz, 3z 2 - r 2 }, {yz, x 2 - y 2 }, {zx, xy}, and ((7 ^ f3)) to 



the pairs {7, (3} — {yz, xy}, {zx, 3z 2 — r 2 }, {?. 



y 2 }- 



e 7 are the on-site energies of the d orbitals. Due to the de- 
generacy of yz and zx, e yz — e zx . From the orbital sym- 
metry it follows £f = t y for 7 = xy, 3z 2 ~ r 2 ,x 2 — y 2 . 
Second nearest neighbor hopping parameters i 7i/ 3 where 
7 = xz,yz and (3 = xy,3>z 2 — r 2 ,x 2 — y 2 change sign 



when 7 and (3 orbitals are exchanged. In any other case 
t-y,p = i/3,7 and i° ^ = , with a = x,y. Other equali- 
ties brought by the symmetry are: 
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The complex sign structure of the hopping terms in- 
cluded in the s x and s y factors arises from changes of 
sign in the orbital wave functions. The factor (— l) m +™ 
in the terms which mix yz, zx with xy,3z 2 — r 2 ,x 2 — y re- 
flect the doubling of the unit cell due to the checkerboard 
alternance of the arsenic atoms displaced up and down 
from the center of the Fe-square plaquettes. These terms 
vanish when the arsenic atoms are in the Fe-planes (see 
Fig. [2] and Appendix I). Due to the enlargement of the 
unit cell, in the reduced Brillouin zone — J < k' x , k' y < -| , 
the Hamiltonian is a 10 x 10 matrix. As discussed in 
the context of three and four band model a 35 i 40 and in 
Appendix II, it is possible to work in an unfolded Bril- 
louin zone — 7r < k Xl k y < tt where k = k' for orbitals 
yz and zx and k = k' + Q in the case of xy,3z 2 — r 2 , 
and x 2 — y 2 ■ In this unfolded Brillouin zone, the system 
is described by a five-band Hamiltonian i?5 X 5(k). The 
relation between the unfolded and the reduced Brillouin 
zones is displayed in Fig. [TJ 

In previous five-band HamiltonianS ) 36 ' 41 the hopping 
amplitudes t% p, p and t 7il g were determined from a fit- 
ting to DFT bands. In contrast, here they are calculated 
within the Slater-Koster framework 3 ^. This method had 
been used before in two and three-band models for iron 
pnictides i 33 ' 42 ' 43 It involves a small number of fitting pa- 
rameters as all the hopping terms depend on a few dis- 
posable constants, the Fe-As and Fe-Fe orbital overlap in- 
tegrals. The final expressions for the hopping amplitudes 
are given in Appendix I. Indirect hopping via As induces 
a dependence of the hopping amplitudes on the angle a 
formed by Fe-As bonds and the Fe plane. This depen- 
dence is shown in Fig. [2j In the range of experimental 
interest of a (29° — 38°) a strong variation of the hop- 
pings is seen indicating important implications for any 
proposed model -either in the wea k 14 ' 22 ' 23 ' 44 or strong 
coupling limi t 33 i 45 ' 46 i 47 to describe these compounds. 

In Fig. [3 and in the rest of the paper the values of 
the overlap integrals and crystal field parameters have 
been chosen to reproduce the main features of the band 
structure of LaFeAsO when a equals the angle measured 
experimentally 4 ^ in this compound a LaFeAs ° — 33.2°. 
While the expressions for the hopping amplitudes given 
in Appendix I include both Fe-Fe overlap integrals up to 
second nearest neighbors, we have found that the band- 
structure is well described including only Fe-As and Fe-Fe 
overlap to nearest neighbors and neglecting Fe-Fe direct 
overlap to second nearest neighbors. With such a choice, 
hopping between Fe atoms to next nearest neighbors is 
completely mediated by As, while hopping to first near- 
est neighbors has contributions from both direct hopping 
between Fe atoms and indirect hopping via As. 

The same values for the overlap integrals and crystal 
field splittings are used in all the figures through the pa- 
per. In section IV we analyze the effect of changing a 
on the band structure, using in the analysis values of a 
which have been found experimentally in several pnic- 
tides. However, we caution on the application of the re- 
sults obtained here on the angle dependence to compare 
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FIG. 2: Dependence of the hopping amplitudes on a. Ex- 
perimental values of a are between 29 and 38 degrees. Top 
(bottom) graphs: first (second) nearest neighbor hopping am- 
plitudes corresponding to pdir=-0. 5, (dda)i = —0.6, (ddn)i = 
0.48 and (ddS)i — —0.1. Direct Fe-hopping between second 
nearest neighbors via (dda)2, (ddn)2 and (ddS)2 is neglected. 
All the energies are in units of (pda) 2 /\ed — e p \ except pda 
and pdir which are in units of pda. Here e p and e<j are the 
on-site energies of the pnictogen-p and the Fe-d orbital (see 
Appendix I). The same fitting parameters and energy units 
are used through all the text. 



different compounds. Substitution of arsenic atoms by P 
or a change in the lattice constant could modify to some 
extent the values of the integral overlaps. 



III. BAND STRUCTURE FOR a haFeAa ° 

Fig. [3] shows the band structure obtained from (1) for 
a = 33.2°, the experimental Fe-As angle in LaFeAsO, 
and the overlap integrals given in Fig. [5J For the crys- 
tal field splitting we take e xy — 0.02, e zx>yz = 0, which 
defines the zero of energy, 63 Z 2_ r 2 = —0.55 and e x 2_ y 2 = 
—0.6, in units of (pda) 2 /|e<j — e p \. The order of the on- 
site energies taken here has been discussed extensively in 
the literature and the values chosen are similar to those 
used previously by other author s 24 ' 36 ' 38 ' 40 ' 49 . The Fermi 
level corresponds to filling the bands with six electrons 
(including spins), as found in undoped pnictides or com- 
pensated FeAs layers. All the figures are shown in the Fe 
or unfolded Brillouin zone. Bands in Fig.[3]have a strong 
resemblance with those obtained from LDA calculations, 
once they are represented in the unfolded Brillouin zone. 
Pockets at the Fermi level include: two hole pockets at 
r = (0,0), a hole pocket in M = (±7r,±7r) and electron 
pockets in X = (±7r,0) and Y = (0,±7r). The resulting 
Fermi surface is plotted in Fig. [5] The two-hole pock- 
ets in r originate in two hole-bands degenerate at the 
top due to the degeneracy of zx and yz orbitals. The 
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(0,0) (n,n) (k,0) (0,ii) (0,0) (n,n) (ji,0) (Q,%) (0,0) 

r m x y r m x y r 



FIG. 3: Band structure in the unfolded Brillouin zone ob- 
tained from the tight-binding Hamiltonian (1) with the hop- 
ping amplitudes computed within the Slater-Koster frame- 
work as described in Appendix I. Values used for the overlap 
integrals are given in Fig. [2] a LaFeAs ° = 33.2°, as found ex- 
perimentally in LaFeAsO. For onsite-energy values see text. 
From (a) to (d) the width of each band-line is proportional to 
its zx, xy, 3z 2 — r 2 and x 2 — y 2 weight. 



so-called Dirac poin t 10 ' 25 which results from the crossing 
of zx and xy derived bands is located close to the Fermi 
level in the vicinity of (0,7r). From the expression of the 
Hamiltonian in Appendix II, it can be seen that these 
two orbitals do not mix in the (0, 0) — (ir, 0) direction. 

The agreement between our results and LDA calcula- 
tions also extends to their orbital character. From top to 
bottom in Fig. [3] we show the energy bands weighted by 
their zx, xy, 3z 2 — r 2 and x 2 — y 2 orbital content, which 
can be compared with the results by Boeri et al£± The 
yz weight is equivalent to the zx weight if X and Y axis 
are interchanged. The two hole-bands in T which cross 
the Fermi level have mostly zx and yz character, while 
some x 2 — y 2 weight can also be appreciated. The elec- 
tron pockets at (±7r, 0)/(0, ±7t) arise from yz/ zx and xy 
orbitals. 3z 2 — r 2 contributes mostly to bands below the 
Fermi level. The orbital content of the Fermi pockets is 
better seen in Fig. [6] 

Around (±7r, ±7r) a 3z 2 — r 2 hole-band nearly touches 
the Fermi level, without crossing it. This band is close 
to another xy hole-band which produces the pockets at 
(±7r, ±7r). All these features are also present in the LDA 
bands. The pocket at (ir, it) has been a matter of discus- 
sion in the literature. Initially-^, it was proposed that 
there was a small three-dimensional pocket with 3z 2 — r 2 
character which, in the reduced Brillouin zone appeared 
at T'. This conclusion was reached using the relaxed 
lattice structure and not the experimental one. It was 
later shown that using the experimental lattice parame- 
ters, in particular the experimental Fe-As angle, the po- 
sition of the top of the xy and 3z 2 — r 2 hole-bands at 
r' in the reduced Brillouin zone switch and the Fermi 




(0,0) (71,71) (7C,0) (0,7C) (0,0) 

FIG. 4: From top to bottom, energy bands corresponding to 
gnashed = 2Q g0 ^ found [q LaFe p ), a rea = 35.3° (regu- 
lar tetrahedron), and a elon9 — 37.2° (elongated tetrahedron). 
The width of the curves is proportional to its xy weight. 

pocket has xy character—. This apparent disagreement 
between the results obtained with the relaxed and the 
experimental lattices originates in a strong dependence 
of the band structure on the Fe-As angl o 13 ' 26 i 27 i 28 . We 
show in Sec. IV that the tight binding proposed here re- 
produces this angle-dependence of the band energies for 
the pocket which appears at (ir, ir). 

The good agreement (shown in Fig. [3]) is not restricted 
to the energies closest to the Fermi level, but it is quite 
generic to all the bands. The correspondence is more 
impressive having in mind that all the hopping ampli- 
tudes are given in terms of just four free parameters. We 
emphasize that the orbital overlaps and crystal field val- 
ues have not been optimized to fit the LDA bands of 
LaFeAsO, but just correspond to the minimum set of pa- 
rameters that reproduce the qualitative features of the 
band structure using the expected orbital energy split- 
ting. 

As discussed above the value of the overlap integrals 
is expected to depend to some extent on the lattice pa- 
rameters and atomic radii. We have found that the band 
structure is sensitive against small changes in the fitting 
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FIG. 5: Fermi surface corresponding to a s «« as ' le < i _ 2g gD ^ 
found in LaFePO), a ^ F <=- 4s0 = 33.2 (as found in LaFeAsO), 
a reg = 35.3° (regular tetrahedron) and a elong = 37.2° (elon- 
gated tetrahedron) with the same fitting parameters as in 
Fig. 3. 



parameters. Close to the Fermi level the largest varia- 
tions appear in the relative position of the top of the xy 
and 3z 2 — r 2 hole bands in M between them and with re- 
spect to those in T, and the energies of xy and zx bands 
in Y. This behavior might be an indication of the exper- 
imentally found strong sensitivity of these compounds to 
modifications in structural parameters. 



IV. FE-AS ANGLE DEPENDENCE 

We now focus on the changes in the band structure pro- 
duced by a modification of the angle a. We assume that 
all dependence enters via the hopping amplitudes. The 
crystal field splitting of the Fe d orbitals results from both 
the As and Fe environment of each Fe atom. Modifica- 
tions in a change the electrostatic environment produced 
by arsenic atoms but not the one due to Fe atoms. We as- 
sume that in the range of a values of interest the change 
in the crystal field parameters is small and we neglect 
the dependence of the onsite energies on a. As discussed 
in Appendix II to analyze the effect of a possible change 
of crystal field with a is straightforward. In Fig. [4] we 
plot the energy bands corresponding to a s i uashed = 29.9° 
(squashed tetrahedron), a reg = 35.3° (regular tetrahe- 
dron) and a elon9 = 37.2° ( elongated tetrahedron) and 
the fitting parameters used in Fig. 3. The bandline width 
is proportional to the weight of the xy orbital. The first 
two values of a used have been found in LaFePO, and 
in BaFe2As2 at optimal doping respectively. The FeAs4 
tetrahedron is slightly elongated in CaFeAsF. 

As evident in Fig. 01 even small modifications of the Fe- 



As angle have an impact on the band structure. Around 
the Fermi level Ep, the most clear change appears close 
to (tt, tt). There are two hole bands with maximum close 
to Ep, with mostly xy or 3z 2 — r 2 character. For a s i uashed 



the 32 



is higher in energy and crosses the Fermi level. 



However the relative position of the two bands changes as 
a increases and the hole pocket around (tt, tt) has xy char- 
acter in the other two cases considered here. The energy 
difference between both bands increases as the tetrahe- 
dron is elongated. This behavior was first obtained from 
LDA calculations. In our tight-binding model it is easy 
to understand the origin of this shift. At (tt,tt) the en- 
ergy of xy and 3z 2 — r 2 orbitals is Ep (tt,tt) — 4t§ + tp for 



(3 = xy, 3z 



The upwards shift in E xy (TT,Tr) with 



increasing a is due to the increase of xy second nearest 
neighbors t xy , while the first nearest neighbors t x xy re- 
mains almost constant (see Fig. [5]). On the other hand 
both t^ z 2_ r 2 and t 3z 2_ r 2 decrease when a increases. 

The dependence of the band structure on the Fe-As 
angle is also seen at T. The gap between the top of the 
hole bands which cross Ep and the xy band at higher 
energies is strongly reduced with decreasing a, in part 
due to a change in t xyjXy (see Appendix II). This gap re- 
duction is, however, not only due to a decrease in energy 
of the xy band. The top of the yz, zx hole bands shifts 
upwards as a is reduced. On the other hand, elongation 
of the tetrahedron can lead to the disappearance of the 
hole pockets at T. 

As a decreases, a transfer of xy orbital weight from 
the third to the second band can be appreciated in the 
(tt, 0) — ► (0, tt) direction. This is accompanied by a shift 
of the Dirac point towards (0, tt). Other changes in the 
band structure with a are discussed in Appendix [Bl 

Somehow weaker, but still observable is the change 
in shape of the electron and hole pockets at T and X 
(Y). This feature is better observed in Fig. [5] The 
electron pockets at X and Y are more elongated to- 
wards r as a is reduced. The shape of hole pockets 
is qualitatively modified as a increases. For the small- 
est angle, a squashed = 29.9°, the hole Fermi pockets 
at r resemble two ellipses centered at T with axis di- 
rected along X and Y directions. With increasing a, 
For both a sc i uashed = 29.9° and a LaFeAs ° = 33.2° the 
two Fermi sheets are very close to each other and would 
be hardly distinguishable in ARPES or quantum oscil- 
lation experiments. With increasing a, for a value cor- 
responding to a regular tetrahedron we find two concen- 
tric pockets. Finally, when the tetrahedron is elongated 
the inner hole has a square-like shape while the outer 
one has a flower-like shape. Similar Fermi surfaces have 
been found in an ab-initio study of the effect of pres- 
sure in the 122 family^. Both circular-like and square- 
like hole pockets at T have been reported from ARPES 
measurements^iii^^^i^^^ii^.. We emphasize that 
we are working in the unfolded Brillouin zone. ARPES 
experiments sample the folded Brillouin zone where the 
pocket that we found at M would be also expected at T. 
Its relative size, compared to the other two hole pockets 
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FIG. 6: From left to right: Orbital content of the Fermi surface corresponding to orbitals zx, xy, 3z 2 — r 2 and x 2 — y 2 . 
From top to bottom, each of the figures is plotted for a aquashed = 29.9° (as found in LaFePO), a LaFeAs ° = 33. 2° (as found 
in LaFeAsO), a reB = 35.3° (regular tetrahedron) and a elon9 = 37,2° (elongated tetrahedron) and the same fitting parameters 
used in Fig. 3. 



in r, will depend on the Fe-As angle. 

In our model, the exact shape of the hole Fermi pock- 
ets found for a given angle can depend on the exact fit- 
ting parameters used. However, the change of the Fermi 
pockets shape with a is a robust feature. 

Recently it has been proposed that the anisotropic or- 
bital makeup of the states on the Fermi surface is crucial 
to determine the superconducting and magnetic proper- 
ties since it controls the anisotropy of the interband pair 
scattering. Fig. [6] shows that this orbital makeup is also 
sensitive to changes in a. Such orbital dependence on a 
will influence the value and anisotropy of the pair scat- 
tering potential. The most dramatic example of such 
sensitivity is the change of the hole pocket at M from 



3z 2 — r 2 , for a = 29.9° to xy, character for larger a 
discussed above. The zx, yz weight in the hole pock- 
ets in T also reverses. For a s i uashed ; the zx weight is 
larger around {hp ,0) in the inner pocket and (0,fci?) in 
the outer pocket. On the contrary, for the regular tetra- 
hedron a reg the situation is the opposite: the zx weight 
is larger around (0, kp) in the inner pocket and (hp, 0) in 
the outer pocket. Note that, in the reference frame that 
we use, zx and yz orbitals lie in the plane of the Fe-Fe 
bonds and are not directed towards the diagonals. Other 
a-dependent effects seen in Fig. [S] include smaller x 2 — y 2 
weight in the hole pockets in T for larger a and changes 
in the 3z 2 — r 2 content of the hole pockets. 



The importance of nesting between electron and 



hole pockets has been emphasized in weak coupling 
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model o ' 22 i 23 ' 36 i which place interband scattering at 
the origin of the magnetic and superconducting prop- 
erties of iron pnictides. Due to the change in shape of 
electron and hole bands with a the amount of nesting will 
be sensitive to changes in the Fe-As angle. According to 
Fig. 5, for the fitting parameter used, the best nesting 
conditions are found between the inner hole pocket in T 
and the electron pockets in X and Y for the regular tetra- 
hedron case. More recently, it has been argued that the 
scattering strength is not simply a consequence of nest- 
ing, but it reflects the orbital weight structure factors. 
The effective pairing strength is larger between fermions 
which belong to the same orbital. In Fig. 6 it can be ap- 
preciated that this nesting is intraorbital, between seg- 
ments of the inner hole pocket and those of the electron 
pocket around Y with d zx character (and around X with 
d yz character, not shown). Interestingly, the same result 
was obtained in ab-initio studies of the evolution of the 
Fermi surface of BaFe 2 As 2 under pressure. 

V. CONCLUSIONS 

In conclusion, we have developed a five-orbital tight- 
binding model to describe FeAs layers in iron pnictides 
with hopping amplitudes calculated within the Slatcr- 
Koster framework. This method to determine the hop- 
ping amplitudes allows to analyze the dependence of the 
band structure on the Fe-As angle a. A good descrip- 
tion of the bands, including its orbital content, can be 
obtained using only four fitting constants to parametrize 
all the hopping amplitudes which compare well with LDA 
bands. 

The flexibility to study changes in the lattice and the 
small number of fitting parameters make this model a 
good starting point to which interactions can be added in 
order to study the magnetic and superconducting proper- 
ties. We have shown that changes in iron-pnictogen angle 
a induce changes in the shape of the Fermi surface and in 
its orbital makeup. In particular, in agreement with LDA 
calculations the hole pocket in (it, ir) (T 1 in the reduced 
Brillouin zone) has 3z 2 — r 2 character for a s i uashed anc l 
xy character for a LaFeAs ° . i n our tight-binding model 
these changes can be understood in terms of the evolu- 



tion of the hopping parameters with a. This sensitivity 
extends to the nesting properties of the Fermi surface. 
In a weak coupling scenario changes of the shape, nest- 
ing and orbital content of the Fermi surface with a could 
be at the origin of the different superconducting order 
parameters, critical temperature and magnetic proper- 
ties found in different iron pnictides. Within the strong 
couplin g 33 i 45 i 46 i 47 point of view the superexchange inter- 
actions will also be affected by changes in a via the hop- 
ping amplitudes. 
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MAT2006-03741 and Ramon y Cajal contracts, and from 
Consejeria de Education de la Comunidad Autonoma 
de Madrid and CSIC through Grants No. CCG07- 
CSIC/ESP-2323 and CCG08-CSIC/ESP-3518. 



APPENDIX A: HOPPING AMPLITUDES IN 
SLATER-KOSTER FRAMEWORK 

In this appendix we give the expressions for the 
hopping amplitudes calculated within the Slater-Koster 
formalism 39 . Both Fe-Fe direct hopping, as well as hop- 
ping via As are included in the expressions below. Fe- 
Fe direct hopping is described via first (dda)i, (ddn)i 
and (ddd)i and second (dda)2, (ddir)2 and (ddS)^ near- 
est neighbors overlap integrals between d-orbitals. Fe-As 
hopping amplitudes are restricted to first nearest neigh- 
bors and involve orbital overlap integral pda and pdir be- 
tween As-p and Fe-d orbitals, which mediate both first 
and second nearest neighbors hopping between Fe atoms. 
To compute hopping via arsenic atoms to second order in 
perturbation theory, we neglect the difference among the 
onsite energies of the d orbitals and that among the onsite 
energies of the p orbitals and take them equal to and e p 
respectively. It is only in the expression for the indirect 
hopping amplitudes that the difference between the on- 
site energies of the d-Fe orbitals has been neglected. The 
values e a are explicitly included and taken into account 
in the tight-binding expression fT]). The resulting finite 
hopping amplitudes are, between first nearest neighbors: 
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xy,xy 



yz,yz 



yz,yz 



— pda 2 — 2pdn 2 + 2V3pdapdir ) cos 4 a sin 2 a + (ddn)i 
e P - e<i | V 2 / 

= ( -pda 2 sin 2 a + Vipdapdir cos 2 a) sin 2 (2a) 

|e P - £d| L V 4 ' / 

+ pdir 2 (cos 2 a + 2 sin 2 a [l - cos 2 a(3 + cos(2a))] )] + (dd<5)i 



1 



f 



32 2 -r 2 ,3z 2 -r 2 



(dd7r)i 
1 

l e P - £ d| 



pder 2 + V3pdapdir I sin 2 (2a) sin 2 a — pd-n 2 (1 — 3 sin 2 a + sin 2 (2a) sin 2 a) 



(Al) 



(A2) 



(A3) 



11 \ 3 

pda 2 sin 2 a ( - cos 4 a — - sin 2 (2a) + 2 sin 4 a J + - pdir 2 cos 2 a sin 2 (2a) 



+ v / 3pc?crpd7rsin 2 (2a) ^— i cos 2 a + sin 2 a 



1 3 
-(dder)i + -(dd5)i 



to 9 9 9 

x z —y z ,x z —y z 



v xy,yz 



3 1 

= -(dd<r)i + -(dd£)i 



1 



t 



yz,3z 2 — r 2 



8V2 
V3 



;«icr 2 sin 2 (2a) + pdn 2 f^l — — sin 2 (2a)^ + pda pdir sin 2 (2 > > ) 



(A4) 
(A5) 
sin(2a()A6) 



^ pder 2 sin 2 a (1 — 3cos(2a)) + pdir 2 y ^ ^— ^ + cos(2a) + i cos(4a) 



+ pda pdir cos 2 (a) ( V2 — — ^= cos(2a) 



sin(2a) 



• a; 

v yz,x 2 —y 2 



1 



3z 2 —r 2 .x 2 —y 2 



1 

l e P - £ d| 



— ^-pdir 2 (l — 2 cos 2 a) — pda pdir cos 2 a 



sin(2a) 



pd7r 2 sin 2 (2a) + pda pdir cos 2 a (l — 3 sin 2 a) 



(A7) 
(A8) 



y^-(dda) 1 + ^j-(ddS) 1 (A9) 



and between second nearest neighbors: 
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^xy,xy 



3 V^3 
--pda 2 cos 2 acos(2a) + pdn 2 cos(2a) sin 2 a — pdapdn sin 2 (2a) 



cos 2 a 



3 1 
+ -{dda) 2 + ^(ddS) 2 



1 



3 "\/3 i X ^ 5 1 

-pda 2 —pdapdn cos(2a) sin 2 (2a) + -pdn 2 II — — cos(2a) — — cos(6a) 

8 2/ 4 V 2 2 



~(eZd7r) 2 + -{ddS) 2 



(A10) 



(All) 



^3^ 2 — r 2 ,3z 2 —r 2 



V3 



»15 \ 3 

pier 2 ( — — cos 6 a + — cos 4 a sin 2 a — 2 cos 2 a sin 4 a + sin 6 a J + -pdir 2 cos(2a) sin 2 (2a) 



pderp<i7rsin 2 (2a) (3 sin 2 a — 1)) 



1 3 
+ -{dda) 2 + -^(ddS) 2 



^x 2 — y 2 ,x 2 — y 2 



l 



1 1 



|e P - e d \ 4a/2 



(A12) 

pcfrr 2 cos 2 a + (dd7r) 2 (A13) 
Zpda 2 cos 2 a cos(2a) + pdn 2 (1 + cos(4a)) - VSpdapdn (cos(2a) + cos(4a)) sin(2a) (A14) 



t 



xy,Sz 2 —r 2 



a/3 

^-pda 2 cos 2 a ( ^ — cos(2a) + ^ cos(4a) ) -r-pdn 2 cos(2a) sin 2 (2a) 



-pdapd-K (1 + 3cos(2a)) sin 2 (2a) 



^j-(dda) 2 + ^j-(dd5) 2 



(A15) 



e P - £d 



3 1 ( 1 1 \ V^3 

- pda 2 cos(2a) sin 2 (2a) pd-zr 2 I lH — cos (2 a) H — cos(6a) I pdapdn cos(2a) sin 2 (2a) 

8 4 V 2 2/2 



+ ^(ddw) 2 - ^(ddS) 2 



1 



2 _„2 = 



tyz,x 2 —y 2 



\e p - e d \ 



— pda 2 (3- 2cos(2a) + 3cos(4a)) + — pdn 2 cos 2 (2a) 
16 ° 



a/3 
2 



+ — pdapdn (cos(2a) — 3 cos(4a)) 



sin(2a) 



, sin(2a) 
2a/2 



(A16) 

(A17) 
(A18) 



Any other hopping amplitude not listed here is zero or is 
related by symmetry to these ones, as discussed in Sec. 
II. The overlap integrals are treated as fitting parame- 
ters. Due to the shorter distance between the atoms the 
largest contribution to hopping is expected to come from 
the first nearest neighbors Fe-As and Fe-Fe overlap inte- 
grals. As shown in the text it is possible to reproduce the 
most important features of the band structure, including 
its orbital content and dependence on a, neglecting all 
the contributions beyond these ones: pda, pdir, (dda)i, 
(ddn)i and (ddS)i. From the above fitting we see that 
while the inclusion of direct hopping between Fe near- 
est neighbors is crucial to reproduce the band structure, 
direct hopping to second Fe neighbors can be neglected. 
Giving all the energies, including d-d orbital overlap inte- 



grals, in units of (pda) 2 /\ed — £p\, and except p-d overlap 
integrals, which are given in units of pda, the determina- 
tion of the hopping amplitudes reduces to the computa- 
tion of just four fitting parameters pdn, (dda)\, (ddn)\ 
and (dd5)i. 

The formalism used here allows to study how the band 
structure depends on changes in the lattice. Indirect 
hopping between Fe-atoms via As induces a dependence 
of the hopping amplitudes in the angle a formed be- 
tween the Fe-As bonds and the Fe-Fe plane. This de- 
pendence is plotted in Fig. [5] All hopping terms show 
angle-dependence with the only exception of t x x i_ y i x i_ y i 
which does not have any indirect contribution. On the 
other hand the amplitudes which couple orbitals yz, zx 
with xy,3z 2 — r 2 ,x 2 — y 2 vanish for a = when the 
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arsenic atoms are in the Fe-plane. 



dk'+Q;yz,(T, rfk'+Q;za;,o-, ^k' ;xy ,a ; ^k';3z 2 - 



r 2 ,<x: 'x 2 —y 2 



the 10 x 10 Hamiltonian can be written as block diagonal: 



APPENDIX B: 



HAMILTONIAN IN MOMENTUM 
SPACE 



H 5x5 (k') 

iJ 5x5 (k' + Q) 



In the ten-orbital reduced Brillouin zone 
— § < k'xik'y < f the Hamiltonian is a 10 x 10 
matrix. Only the terms which mix orbitals zx,yz 
with xy,3z 2 — r 2 ,x 2 — y 2 feel the unit-cell doubling 
and couple states with momentum k' with other with 
momentum k' + Q with Q = (tt, tt). Choosing the 
orbital basis dk' ; a,<T, dk'+Q;a,<r in a convenient order: 



k';yz,ai dk' lZ x,cr, +Q;xy ,cr , CEk'+Q;3z 2 -r 2 ,eri a k' +Q;x 2 ~y 2 ,ai 



with 



H 5x5 (k') = 



H 2x2 (k') 
H 3x2 (k') 



H 2x3 (k') 
^3x 3 (k' + Q) 



- (fx - e 7 )I 



Here 1 is the unit matrix and 7 the orbital index in the 
order given above. The subindices in the matrix name 
serve to label the matrices and indicate their dimension. 



#2x 2 (k') 



2tyz,yz C0S ky 



-4L 



'yz,yz COS k x COS ky 

^yz,zx sin k x sin ky 



^yz,yz C0S ^x 



^ V zx,zx C0S ky 



-& V z,zx sin k' x smk!y 

+ 2t zxzx COS k x + tt zx ^ zx 



cos k' x cos k' y 



I 2t % y ,xy( COsk 'x+ COsk 'y) 
~\~^t>xy ^xy COS k x COS ky 



H 3x3 (k') 



-tt 



xy,3z 2 —r 



_ r 2 sin k'. sin k{, 



V 



-tt xy3z 2_ r 2 sin k' x sin k' y 



_ r2 (cos k' x + cos k'y) 
4i 3z 2_ r 2 3z 2_ r 2 cos k' x cos k' 



2t x 

3z 2 —r 2 ,3z 



2t 



x 2 —y 2 ,3z 2 —r'' 



(cos k' x — COS ky) 



(cos fc^, — COS k'y) 



2t x 

x 2 —y 2 ,3z 2 —r'- 



2t x x2 _ y2 x2 _ y2 (cos k' x + cos k'y) 
+tt x 2_ y 2 x 2_ y 2 cos k' x cos k'y / 



and H 3X 2 = [^2x3] T with H 2x3 (k') equal to 



2i sin k'y (t v yzxy - 2 t „ , , x „ c o s fc ^ ) 2isinA;U^ ;Zi 3 Z 2_ I ,2 -2f J/Zi3z 2_ r .2 cos fcQ 2i sinfc^ \t x yz>x z_ y ?. -^yz,x 2 ~y 2 cos fc^) 



2z sin £4 (t zx xy 



^t-zx^xy COS <i/^. 



y \ ^zx,3z 2 — r 2 



-2t 



zx,3z 2 



cos k') 2i sinfc' ( t 



il 



y \ zx,x 2 —y 2 



^zx,x 2 — y 2 COS k x 



In the above mentioned basis order, the block diagonal 
form given provides a natural wa y 35 i 40 to unfold the Bril- 
louin zone (see Fig. [T]), that is, to define k = k' for or- 
bitals yz and zx and k = k' + Q for xy, 3z 2 — r 2 and 
x 2 — y 2 . In this extended Brillouin zone — tt < k x , k y < tt 
the Hamiltonian is given by i?5 X 5(k). All the figures and 
expressions in the main text are given in the unfolded 
/c-space. 

One advantage of the present tight binding model is 
the possibility to understand many features of the band 
structure. In particular exactly at T and M none of the 
five bands show orbital mixing and simple expressions 
follow for their energies: 

£ 7 = ±2t» 7 ± 2t x lrl + 4t 7 , 7 + e 7 - fi (Bl) 

with 7 = yz,zx. Plus (minus) sign applies at T (M). 

E. ( = T 4i^ 7 + tt~ ln + e 7 - n (B2) 

for 7 = xy, 3z 2 — r 2 and x 2 — y 2 . Minus (plus) sign applies 
at r (M). From (|B1[) and ||3J| the degeneracy of yz and 



zx bands at T and M follows. This two-fold degeneracy 
is clearly seen at the top of the hole bands which cross 
the Fermi level in T and in the two highest in energy 
bands at M. As discussed in the text the dependence 
of E xy (M) on a originates in the sensitivity of t xy ,xy to 
changes in the angle. The same dependence is present in 
E xy (T) which shifts with a in the same way as E xy (M) 
does, namely, decreases as the angle is squashed, keep- 
ing E xy (T) — E xy (M) = —8t xx almost unchanged. On 
the contrary, due to the combined dependence of first 
and second nearest neighbors in a both E 3z 2_ r 2(M) and 
E 3z 2_ r 2 (T)—E 3z 2_ r 2 (M) decrease when the angle is elon- 
gated, the latter becoming eventually negative. Depen- 
dence of the hopping parameters on a are plotted in 
Fig. 2. 

At X and Y only 3z 2 — r 2 and x 2 — y 2 mix, the energies 
of the other orbitals can be expressed in the simple form: 

E 1 (X, Y) - ±2t» 7 T 2t x y - 4f 7 , 7 + e 7 - M (B3) 

for yz and zx. First sign applies for X and second for Y . 
The degeneracy found in T and M is broken, but their 
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energies are related by symmetry E yz>yz {X) — E ZX . ZX (Y) 
and E yzyz (Y) = E ZXZX (X). Due to the sign which 
precedes t xy>xy in the expression for the energy corre- 
sponding to the xy orbital in X and Y, E xy (X,Y) = 
—4&xy,xy + e xy — fi>, the shift of E xy (X, Y) with a is op- 
possite to that found at T and M. 

In the present paper we have neglected the dependence 
of the crystal field splitting on the angle a but, due to the 
simplicity of the expressions for the energy at the sym- 
metry points, guessing its effect in the bandstructure is 
straightforward. In particular, a possible change in crys- 
tal field of xy with a would shift E xy in the same amount 
in r, M, X and Y, contrary to the effect produced by 
the angle-dependent hopping parameters. 

Another interesting feature regards the mixing be- 
tween orbitals along the high symmetry lines TY, TX, 



MX and MY. Along FY and MX xy/yz bands cross 
zx/x — y 2 /3z 2 — r 2 bands without hybridization result- 
ing in Dirac points. These Dirac points can be observed 
in Fig. 3 and Fig. 4. Along MX there are crossings be- 
tween the two upmost bands and between the second 
and third bands (the later crossing being only present 
for a s i uashed i n Fig. 4). Along TY such crossings ap- 
pear between the two upmost bands and between the two 
lower ones, the former crossing being absent in the case 
of a s i uashed , Along this direction a Dirac point close to 
the Fermi level is also found at the crossing between zx 
and xy derived bands, as discussed in the main text and 
mentioned previously by other authors i 10 ' 25 The same 
physics appears along TX and MY with the interchange 
of yz for zx. 
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